Iran. J. Insur. Res., 12(3): 197-212, Summer 2023, Serial #45 


Iranian Journal of Insurance Research 
(WIR) 


® 


Insurance Research 
Center (IRC) 


Homepage: https://ijir.irc.ac.ir/?lang=en 


ORIGINAL RESEARCH PAPER 


Pricing of life insurance products using markovian aging process model 


A. Rostami*, A. HassanZadeh 


Department of Actuarial Science, Faculty of Mathematics, Shahid Beheshti University, Tehran, Iran 


ARTICLE INFO 


ABSTRACT 


Article History: 

Received 12 February 2023 
Revised 05 May 2023 
Accepted 04 June 2023 


Keywords: 

Aging process 
Markov process 
Mortality prediction 
Physiological age 
Pricing 


*Corresponding Author: 
Email: ar_rostami@sbu.ac.ir 
Phone: +9821 88554101 
ORCID: 0000-0002-3754-420X 


DOI: 10.22056/ijir.2023.03.03 


BACKGROUND AND OBJECTIVES: In this research, our main objective is more accurate 
pricing of life insurance products with a new approach of predicting mortality or survival 
rates. Currently, a life table is used to calculate the current value of pensions, insurance 
premiums, etc. Therefore, to increase the accuracy of our calculations, we are looking for 
a mortality prediction model for such calculations. Therefore, in this research, instead of 
static pricing (only using the latest life table), we used life table prediction and dynamically 


rated life insurance products. 


(METHODS: In this research, a new model proposed to predict the probability of human 
mortality (survival) based on the Markov process, a limited state with an absorption state 
(death). This model measured based on the physiological age, because the physiological 
age of each person can be checked based on different laboratory indicators, and finally 
it has led to the results of the individual health index. In addition, the parameters of this 
model are the initial probability vector and the sub-intensity matrix of a Markov chain 
that changes over time. In other words, in this model, according to a possible process 
in the model, the initial probability vector over time selects the possible interval of the 


physiological age equivalent to the chronological age. 


FINDINGS: To show the satisfactory performance of this model, the relevant data set from 
the United States of America was analyzed. The predicted results with the presented model 
are better than Lee Carter’s model. It should be noted that the number of parameters of 
the model introduced in this research is much less compared to the Lee Carter model and 
other mortality or survival prediction models. Based on this model, a closed form for life 


insurance pricing relationships is obtained, which simplifies these calculations for users. 


CONCLUSION: The relationships obtained for pricing were investigated based on two 
products, 5-year term life insurance and also a 5-year term pension. The fitted results for 
the model used in the predictions of the probability of mortality as well as the probability 


of survival and pricing are very satisfactory. 


This is an open access article under the CC BY license (http://creativecommons.org/licenses/by/4.0/). 
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Fig. 1: Proposed Markovian physiological aging process 
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Fig. 2: log of the probability of death between the ages of 30 and 108 and from 1970 to 2014 in the USA 
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Table 1: Estimation of parameters 
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Fig. 3: Mortality index varying with time k(t) 
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Table 2: Time series model for forecasting 
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Fig. 8: Forecast of death rate for different ages from 2015 to 2019 
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Table 3: Comparison of the actual value and the forecast of life expectancy for the years 2015 to 2019 


2017 


Ls 
2019 2018 2017 2016 2015 J 
Year 
: SSH) & deel 
ily . 
50.0 49.8 49.8 49.8 49.7 ALL Ys 955 
Actual : 
Life 
oe expectancy of 
50.4 50.3 50.1 50.0 49.9 iiaciiad a 30-year-old 
Forecast person 
: Sy Fj A eal 
stils ; 
31.7 31.5 31.4 31.4 31.2 aJLs d+ 958 
Actual 
Life 
i, expectancy of 
31.8 31.7 31.6 31.7 31.3 ei a 50-year-old 
Forecast 
person 
. SF) & deel 
15.8 15.6 15,5 15,5 15.4 ee dita¥eas 
Actual 
Life 
er expectancy of 
15.8 15.7 15.6 15.5 15.4 iis a 70-year-old 
Forecast 
person 


Y-Y 


YIV-VVV Ge AREY Glee LP 0 laut AV 0590 dow doling} cole ay pti 


Sete 9 ILE Ogle eye FF Jone 
Table 4: Comparison of goodness of fit and prediction 


. : les 
Se lis le pal, bolas 
or shat 9 poe do 
ee Survival BIC AIC Numberof “2? "~” MSE R2 J 
Mortality Log- likelihood Model 
forecast parameters 
forecast error 
error 
x he ae | Ke 
0.0277 0.0161 36080.2 35771.4 50 -17835.7 0.0209 0.99991 eon ol J 
Our Model 
PIS J J 
0.0632 0.0292 37125.7 35884.3 201 -17741.2 0.0139 0.99994 Lee-carter 
Model 


WS ge CSI Mia Gl! a5q0,@ lbF dan Sd 43 50 | 
Lei pplaio cladan G> b clrolo 8 Gly cul SS & p5¥ 
peels Gel G2 Til thy he ks ae Ble es 
LL 59 g dF dan 98 Glo} 59 road Ogu Ke fold blir. 
3 eSB 9 - May GLI og5 45 plas ye - ods are ole} 
98 NF do ST cel shite Ip So 9 ile} pos don 
Sl 9 99,5 ce Sly ile poe don Geb dlepe oly! O98 
hot my paw (yloy 9 Sho y1,3 Bolles ale pw aisle, 085; JbF aay 
WS ce II) plate CSO Sy Gpotno closlo 8 SQ 
pri ON pF 958 dq} 4 Cpetre CSloyp AF SL; 
OPIPD 9 Opt ce ctrl poe Gyre qs poime adh arilo 
(Poterba, 2001) o9 co Anat Sac poten cod SSL yo 4 
LrSy Slop (bed G55)1 eb iI LE ee aslo! jo Gre 
Ole age Cal edb 8 ne Joe lel 1) Cpotme Gaizem 9 
p35 (52 
owl C (Bolai CaS chlo gu Gisy 99 od 320 Jo 
Cy 0dgdrre AS iS co ole 1) Slob WS a oll cyl 
balgy cal ply seocge Garetd 55 Ly Je 39 99 Ride d 
gp SBlgs pj S ygrds Le Jae oluly 6 SaaS 
15398 Glo} 59 LEX Slop LS (85,1 © 
f e~Sa,e7*(—T,)1ds @) 
i-1 
= A.og = Le, 0A8I— Ty fe 8 — g - —-2) 1 


a 


i 
axl = | e~Sa,eT1ds el 
i-1 


i=1 


= Gem j=1 4; (41 — T,)1{e 6 YE — ¢GI-M)}4 


Y- 


PWS I Joe jl pay le Jao 59 AIC 9 BIC az — Col 5,15 
BIS I Je celsius ze pW 9 MSE aie jo a! cu 
pak obj class egdg0 gel cde af 8b ge Joe I te 
AS Sal Se jo cpl cl le Jae Gen FIT J Jue 50 
eib ce be Jae ply go pte SIT I Joe Grote lke 
dee BIS I Joe jl Ae be Joe 5 ete Oo ale 

ABS, os 


(rw pote 9 Oy T 

Lm jgdS adn 8 idle E55 5B ej CES i! © 
ple] faa Gise Gal go Ceol ppttine Ghulro Gynizred 9 
me 9 S5tce AIL om 9 6 pSgAST clad (Sj)! droleee 590 59 
BaF vs Bh) Slabs pales 

Vgane Gpoime pot Gade 9 Ab Sogo Glacsh y 
Sglrre slay Cle 59 9 op co ples! plas kj cleojl 59 
g ee SV gare CT 85)! sober GL Spee The ol 
Lode Ge Slslrs jo Se Cos) Cpesl le jl C pote 
Lisl |) GAB Gpasiege (cSimtijs colalja ao 9 lndaliden «85,1 
Cel Bola pttie Kyat Gpoine Ke hd 55,1 UTE 
Gyoine poe don claslo 1,3 0 p)o Kim 958 poe Job a a1 > 
dow ABbdewgs clay oled jo 099) delllacd go Eqdge aS 
eek Sl poe don GV pare cole 5 jb oj) og>5 
9 YL claolo 5 Gate Ce pe jLS0 90 CAKES 9 Cul 00 SF 
(Mitchell, 2002).cus) ord Sorren 6X0 le) > 5 (jydo 
SY gare Lhre pc don 9 pocpled don hj poe don 
Sewype 9 Lg S38 Ole} yo GAH Cadre fold a5 alc 
pot dow db 0 Ce je 9 dow Go odd pei hj] polde L 
Oe Oddday C98 Glo) 99 ST 9 0d Codd Ogu SQ Jbj 
Aes iS go Sadly tates Glin! yadete Clee Si ghly 3 
Syed CoE lade So E aaoice ojlel iF don 4 ile} pos 
erred ale pacpls dan lei pals gloolsil cl os 


A 


IN jhe OY pare 6 IF O05 


Siglitie (cleo ys Ey g Gai oly Aggy Gh sgeST bd 95,1 0 Jgue 
Table 5: Present actuarial value of 4x:5| for different ages and interest rates 
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